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ABSTRACT 


An alternate description of the Plasma Opening Switch 
(POS) is presented. In addition, a preliminary study was 
conducted into the ionization of desorbed carbon particles 
from the cathode surface in the switch region. The formation 
of a locally dense plasma layer surrounding the cathode is 
Shown to influence the electric field and the Debye screening 
distance. This effect influences the emission of electrons 


from the cathode as well as providing a mechanism for opening 


the switch. 
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TL. INT ROBUEION 


A. BACKGROUND 

Inductive energy storage systems are capable of storing 
energy at densities up to 100 times that of capacitive 
systems [Ref. 1]. However, inductive storage systems are 
dependent on their opening switch characteristics. Switches 
range in type from mechanical (slow) to exploding fuses 
(faster) to plasmas (fastest). 

fiemwvomeage OouLput provided by an inductor is given by 
V = -L di/dt, where L [H] is the inductance and di/dt [A/s] 
is the changing current. As the switch determines the di/dt 
term, it is clear that faster switch performance can provide 
larger voltage/power output with smaller input current 
sources. 

Currently, plasma switches are being incorporated into 
vacuum inductors to provide the fast switching necessary for 
pulsed-power applications. Of particular interest is the 
Plasma Opening Switch (POS). The POS has been used in 
experiments with the Gamble I [Ref. 1] and the Gamble II 
[Ref. 2] accelerators at the the Naval Research Laboratory 
(NRL) as well as with the Pollux and Kalif pulsed-power 
generators [Ref. 3] at the Kernforschungszentrum in 
Karlsruhe, West Germany. The results of these experiments 


have demonstrated that the POS was capable of conducting very 


large currents before opening in less than 10 ns. Meger 
[Ref. 1] reported voltage and power multiplication of 2 and 4 
times that of experiments without the switch while Ottinger 
reported Similar results from experiments using the Gamble II 
accelerator [Ref. 2]. The POS was renamed the Plasma Erosion 
Opening Switch (PEOS) when Ottinger presented his theoretical 
model of the switch [Ref. 4]. 

The Naval Research Laboratory has been actively using the 
PEOS to provide pulse compression for their experiments 
involving a plasma radiation source (PRS). These 
experiments, sponsored by the Defense Nuclear Agency (DNA), 
were undertaken to develop a source which simulates the 
radiation produced by nuclear explosions. The results of 
these experiments have demonstrated that the performance of 
the PRS was directly affected by the use of the PEOS. The 
PEOS significantly reduced the prepulse and compressed the 
power pulse rise-time. This effectively inhibited the growth 
of instabilities in the plasma radiation source and produced 


a more efficient yield. [Ref. 2] 


B. THEORETICAL MODEL OF THE PEOS 

A theoretical model of the PEOS has been developed by 
Ottinger [Ref. 4] consisting of four phases: conduction 
phase, erosion phase, enhanced erosion phase and magnetic 
insulation phase. A complete description is available in 


Reference 4. 


During the May 1986 DNA Pulsed-Power Review Meeting in 
Las Vegas, Nevada, the Maxwell Laboratory presentation 
reported the results of their latest PEOS experiments on the 
BlackJack 5 generator. The results demonstrated PEOS 
Switching of currents up to 2 mega ampere with prepulse 
Suppression. However, the report also indicated that the 
Switch opened too slowly for pulsed power applications and 
that the magnetic insulation of the cathode surface was 
imeenpleee.weelae Conclusions stated that critical aspects of 
the switch parameters and the physical mechanisms of the PEOS 
model were not well understood. [Ref. 5] 

Schwirzke has asserted that the ion diode modeling of the 
PEOS is not applicable and the physics describing the plasma 
Switch are incomplete. He has proposed an alternate 
description of the physical processes of the switch based on 


reduced field emission of electrons from the cathode surface. 


C. PLASMA SHEATH MODEL 
The proposed model is based on Schwirzke's earlier work 
[Ref. 6] and the recent work of Halbritter [Ref. 7]. These 
experiments demonstrated: 
x Plasmas in contact with metals induced significant 
desorption of gases and other impurities from the 


surface. 


x Electron field emission at contaminated surfaces is 
achieved with reduced electric fields. 


These observed results led Schwirzke to postulate an 


alternate description to the PEOS model. 


Schwirzke has asserted that Ottinger's PEOS model does 
not fully describe the physics G@ee@hemewitch- =schwiezke's 
major points of disagreement are: 


x The electron current is field emitted rather than space- 
charge-limited. 


x The sheath that forms between the bulk plasma and the 
cathode surface adjusts to the applied voltage at the 
cathode rather than being fixed. 

x There is no physical mechanism which causes the PEOS gap 
to increase. 

BD. PURPOSE OF Tape Reo ean? ih 

The purpose of this research is to qualitatively describe 
Schwirzke's opening switch model and perform some preliminary 
testing to explore the influence of the neutral carbon 
ionization on the performance of the switch. Other aspects 
of the model will be treated gualitatively. It is hoped that 
this research will serve as a basis for the continued 


development of this model and that important insight into the 


Switching behavior of the plasma can be achieved. 
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Il. PLASMA SHEATH MODEL 


A. OVERVIEW 
1. Purpose of the Model 

The proposed plasma sheath model is an alternate 
description of the switch model known as the Plasma Erosion 
Opening Switch which has been developed by Ottinger [Ref. 4]. 
The plasma sheath model is based on the reduced field 
emission of electrons from a contaminated cathode surface. 
This model attempts to simplify and gain insight to the 
understanding of the physical behavior of a highly non- 
uniform plasma within a confined switch region. 

22a Opening Switch Description 

The basic configuration of experiments using the 
plasma opening switch (POS) is shown in Figure 1. It 
consists of a generator section, vacuum coaxial inductor 
section, switch section and load section. Plasma is injected 
from plasma guns located off-axis toward the cathode. After 
a delay period, the generator is fired and current begins 
flowing from the outer anode to the inner cathode through the 
plasma switch region. The switch region conducts the current 
flow for some period of time before opening quickly and 
providing a high voltage pulse to the load section. It has 
been shown that the performance of the opening switch is 


highly dependent on the plasma density, the delay time for 


at 


firing the generator and the physical dimensions of the 


Switch region. [Ref. 8] 
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Figure 1. Experimental Setup 


B. PLASMA SHEATH DESCRIPTION 
1. Assumptions 
The major assumptions which apply to this model are: 


* The anode-cathode potential difference is applied over 
the sheath region. The effective resistance of the bulk 
plasma is very small and it is determined by the Spitzer 
resistivity and the switch region dimensions. The 
potential difference required to drive the current 
through the bulk plasma is minimal and can be neglected 
compared to the total anode-cathode potential 
difference. 


x The stainless steel cathode has several monolayers of 
contaminants on the surface. The main source of the 
contaminants is from the carbon plasma guns. After each 
firing of the plasma guns, a certain amount of residue 
settles on the surface of the cathode. For the purposes 
of this research, only the carbon neutrals will be 
considered. 


WZ 


The streaming gun plasma induces desorption of the 
cathode surface contaminants. The desorbed contaminants 
expand, forming a dense layer or cloud of neutral 
particles. These particles are ionized by the emitted 
electrons from the cathode surface. 


esiitetalecurrent tlow Gmethe Swatch region 1s uniform 
and in the radial direction. The field emitted 
electrons, normally emitted at very small sites such as 
whiskers or dirt spots, are considered to be distributed 
over the entire cathode surface in the switching region. 
In addition, the ion flow is considered to be uniform 
with a constant drift velocity toward the cathode. 


The switch dimensions and the density of the bulk plasma 
remain constant throughout the duration of the switching 
period. The density of the bulk plasma is confined 
without losses through any boundaries with the exception 
of the ion current at the cathode surface. However, 
this loss is small and can be neglected so that the 
plasma density is essentially constant. 


The sheath region does not break down due to the 
voltages across the gap. In addition, the plasma 
temperature assumes the energy of the cathode. As the 
cathode energy increases, the sheath distance increases 
by pushing the electrons further away. This imparts 
energy to these electrons which in turn is applied to 
the bulk plasma. 


The cathode is cold so that thermal emission of 
electrons is not significant and can therefore be 
neglected. 


Desorption 


As the injected plasma comes into contact with the 


cathode, the neutral contaminant layers are disrupted and 


these atoms are desorbed from the surface. These neutrals 


expand slowly (approximately 330 m/s) and form a locally 


dense layer of atoms to be ionized. The expansion of the 


neutral particles is considered only during the time period 


that the plasma is in contact with the cathode surface prior 


to the firing of the generator. This is due to the slow 


is 


expansion velocity and very short time frame for current 
conduction. Therefore, the number of neutral particles 
desorbed and the distance traveled from the cathode surface 
is considered constant throughout the switching phase. 

3. Initial Field Emission 

The initial emission of electrons from the cathode 
surface is produced by localized spots with high electric 
field strengths. The electric field, E (V/m) is enhanced by 
a factor & such that the enhanced field is given by F = RE. 
The enhancement factor range is 10 € &® S$ 1000. 

Halbritter's work {Ref. 7] demonstrated that the 
electron field emission current from contaminated surfaces 
was observed at electric field strengths on the order of 
F = 10° V/m epithe contributing factors to field enhancement 
are dust particles, whiskers and other contaminants. 

The electric field is defined as E = - Vo, where 0 [V] 
is the potential field. For the case of cylindrical 


coordinates and a potential field which only varies in the 


radial direction, then E = -d@/dr. Furthermore, if it is 
assumed that dd/dr is constant within a given region, then E 
can be approximated by the potential difference over that 
region. 

The plasma sheath model assumes that the anode-cathode 


potential difference occurs within a sheath region given by 
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AD= 2 (1) 


Wiebe we. |F/m) 1S the permittivity of free space, e [C] is the 
elementary charge, 9. [V] is the cathode potential and n [m>] 
is the electron density of the plasma. This is the Debye 
screening distance for a negatively biased cathode. As 


stated earlier, E ~ 0/A, so that E varies with the plasma 


density and cathode potential as 


1/2 
ne o. 
E= (2) 


€ 
O 





ieee Pomnzat On 

tie stenmizatlon rabewor the neutral particles is 
dependent on the neutral density, ionization cross section 
and the number of incident electrons. This research 
considers only the neutral carbon atoms which are ejected 
PoMeeieesurrace fOr Calculations of the ionization cross 
SeCE LON ; 

As the neutrals become ionized, a locally more dense 
plasma forms near the cathode surface which influences the 
electric field and the sheath width. From Equation 2, it is 
clear that for a constant potential, the electric field 


varies with the density. For increasing density, the 


ies 


electric field strength increases and additional field 
emission from the cathode is possible. 
5. Ion Removal 

The ionization of the neutral particles produces a 
locally more dense plasma near the cathode surface. The 
1onized particles in the sheath, under the influence of an 
enhanced electric field, are accelerated toward the cathode 
surface. The time required for these ions to reach the 
surface is determined by the accelerating field and the 
distance traveled. As the carbon ions are removed from the 
sheath region, they contribute to the total ion current flow 
due to the higher density. In addition, the reduction in the 
local charge density causes the sheath width to increase and 
the electric field to decrease. As the magnetic field 
strength increases, the electrons begin to drift toward the 
load side of the switch before crossing the gap and entering 
the bulk plasma. 

6. Magnetic Insulation 

Magnetic insulation occurs when the magnetic field 
near the cathode is sufficiently large to cause the electron 
Larmor radius to be smaller than the sheath width. As the 
magnetic field experienced by the emitted electrons 
increases, the electrons experience the ExB (E iS assumed to 
be radial and B is in the @9@ direction) drift which causes the 
electrons to drift parallel (z direction) to the cathode 


surface for the length of the switch. The current is 
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switched to the load and the effective switch resistance 
increases. The switch is considered open when the Larmor 


radius becomes less than the sheath width. 


Ey 


Lid. JONITIZAT IONSO RANE Ural PARTICLES 


There are four critical areas of interest in the plasma 
Sheath model which need to be addressed. These are: 


x The ionization of neutral particles which are desorbed 
from the cathode surface. 


* Electric field strength and sheath width response due to 
the increased ionization of neutral particles. 


x Ton removal occurring at the cathode surface. 
x Magnetic insulation which causes the switch to open. 
This section is devoted to the calculations necessary to 

determine the number and density profiles of the desorbed 
carbon atoms, the number and density profile of the newly 
formed plasma and the electric field produced by the pressure 
force associated with a density gradient. The computations 
used are based on fundamental plasma physics relationships. 
No effort has been made to apply statistical methods to this 


problem. 


A. DESORBED =NEULRAES 

The cathode used in the experiments at the Naval Research 
Laboratory and at the Kernforschungszentrum in Karlsruhe is 
made of stainless steel. Most stainless steels consist of a 
small percentage of nickel and chromium (usually below 20%) 
mixed in with the base metal of iron. At temperatures below 


912 degrees Celsius, the crystal lattice structure of iron is 


ls 


body centered cubic (BCC) with a lattice parameter of 
2.866x1077>° meters. The contaminants, which are situated on 
the surface, will occupy the interstitial sites. Therefore, 
the BCC configuration and lattice parameter of iron will be 
used for the calculation of the number of carbon particles 
desorbed from the cathode surface. [{Ref. 8] 

The volume of a layer on the cathode surface within the 


See reg@on 1S given by 
eS 27mR I cr (3) 


where R. [m] is the radius of the cathode, 1 [m] is the length 
of the switch region and dr [{m]} is the depth of the surface 
layer. The value of dr is given by dr=0.5*m*2.866x10 ?°, 

where m is the number of monolayers desorbed. The BCC 

lattice contains two atoms/cell, therefore the number of 


atoms desorbed is given by 


2nR_| m (4) 
= ———___ 
(Dede 


Be Se NBUTRAL DENSITY PROFILE 
The radial expansion of the neutral carbon layer is 


determined by 


else rt V (2) 


ARS, 


where R, [mM] is 9Gne OfE axis IoeaetiGn@er Chieu p fant -mcnunio mae 
[s}] is the delay time from plasma gun firing to generator 
firing, vg [m/s] is the plasma deimitevelocity anda =m olen. 
the neutral expansion velocity. The neutral particle 
expansion distance will be considered constant for the 
duration of the switching action due to the short time frame 
involved and the relatively slow expansion velocity of the 
neutrals. 

Sub-regions are formed by dividing Equation 5 by the 
number of sub-regions desired. Initially, each sub-region is 
assigned an equal number of particles based on the 
calculation of Equations 4 and 5. The profile can be 
adjusted so as to weight the number of neutrals in a given 
sub-region as will be explained in Appendix A. From this 
number profile, a neutral density profile for each sub-region 
will be calculated from 


N 
p= (6) 


fs Z 
m(R-R. ) 
where Np [m*] is the number of particles per sub-region and 
R 


5 (m] and R; [mm] ane Ehe outer anemone det someGin oe Cavon 


sub-region. This neutral density profile is important for 


determining the mean-free-path of the electrons within each 


Sub-region. 


ZO 


C. BBECTRON FLUK 

The number of electrons flowing through the layer of 
desorbed neutrals is determined from the electron current 
density flowing through the cathode surface. For a given 


tomeal Current, I [A], and a given 10n Current, 
i= nev AG (7) 


where A, [m*] is the cathode surface abea, the electron flux 


at the anode is determined by 


[ - I. es 
[= eA. “a 10 (8) 





where A, [m7] is the anode surface area and I has the units 


[number of electrons per square meter per nanosecond]. The 
electron current is extracted from the anode surface while 
the ion current is extracted at the cathode surface. 
Conservation of charge requires that the electron current 
extracted at the anode must equal the electron current at the 
cathode. Therefore, the radially outward electron flux into 


ecylindrical region is given by 


r=Tl Rp (9) 


a 


where R, [m] is the radius of the inner surface that the 


electrons flow through. The number of incident electrons per 


nanosecond into each cylindrical region can be expressed by 





I- I ” 
Be |e (10) 





which is constant for a given current. 


D. ‘DENSITY OF VIONTZEewE ART ICEES 


The mean-free-path of the electrons is determined by 
-1 
A=(n.0) (11) 


where GO [m*] is the electron impact ionization cross section 
for the given target material. The value of O is dependent 
only on the energy of the incident electrons. Table 1 is 
used to determine the electron impact ionization cross 
section for neutral carbon atoms [Ref. 9]. The values used 
range from relatively low energies ~ 20 eV up to ~ 1000 eV. 
It is assumed that the ionization occurs at low voltages so 
only relatively low electron energies will be considered. 
The values of energy are given in electron volts and the 
cross section is given in square meters. 

The ionization of the neutral particle cloud is dependent 
on the electron flux, mean-free-path and the thickness of the 


medium. For a given flux of electrons incident upon a medium 
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of thickness x and mean-free-path, the flux of electrons 


exiting the medium is reduced by the factor exp (-x/iX). The 
difference between the incident and exiting electrons 


represents the ionization that occurs in the material. 


TABLE 1 
Penet RelmeMPACT SE@NIZATIONSCROSS SECTION 


Energy Cross section Energy Cross section 
20 ths QUOT 210 200 1 SON SAO) 
30 oe OZ} 300 eee — 20 
40 2 MO SZ2.6 400 oe Oa 
50 EOE 20 500 SOO =2 0 
60 2 405-20 600 7. OUR 21 
70 Ze 0b - co 700 Caen —2 1 
80 BEN 2.0 800 one OB—Z 1 
90 Ze tee) 900 Ss). SO 2 al 
100 Zoi 2 © 1000 Se Oe al 


Tie number of neutral carbon atoms that are ionized can 


be expressed by 
N.= nN | 1 - exp (-x/2) | (12) 


It follows that the density of the desorbed carbon atoms that 
become ionized is given by 
N. 


n.= ee (13) 
ml (R2- R?) 


At the end of each nanosecond time increment, the number 


of particles that become lonized in each region is subtracted 


from the original number. The new densities, mean-free- 


Vee 


paths, etc. are then recalculated for the next nanosecond 


increment. 


BE. DENSITY RATIO AND PRESSURE GRADIENT 
The ratio of the ionized carbon particles to the plasma 


density is ogven by 


n 
D=— (14) 
n 
where nm, [m7>] =n, + n, is the total ionized density. The 


electric field at the cathode and the sheath width are 
influenced by increased values of D. The electric field 
experienced at the cathode can be modified by multiplying 


\/@ which demonstrates the influence of the 


EQuation. 2 .by 9) 
ionization on the elecurmicerield: 

In addition to the influence of D on the electric field 
experienced at the cathode, an increasing D produces a 
pressure force within the plasma. Treating the plasma as an 
ideal gas, the pressure can be expressed as p = nkT, where k 
[J/K] is the Boltzmann factor and T [K] is the electron 
temperature of the plasma. Considering only the pressure 


gradient force, the electric field produced by the gradient 


can be expressed by 


=n (i>) 
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where ng [m7>] is a modified density term. This term is the 


average density between adjacent regions and is given by 





(16) 


if the pressure gradient force is sufficiently strong, the 
motion of the ions and electrons in this region will be 


arrecred. 
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iV.  BeSUinisS 


The results of this study are presented in two parts. 
First, the 2tonization occurring @meene desorbed carbon cloud 
is evaluated. Second, the effect of the pressure gradient on 
the electric field due to the ionization is considered. The 
Switch parameters and the input current used in this study 
are taken from the Gamble I experiments reported in Reference 
10. The effect of varying the switch parameters on the 
ionization and resultant electric field strength is also 
evaluated. 

Typical switch parameters given in References 1, 2, 4 
and 10 are shown in Table 2. 


TABLE 2 
SWITCH PARAMETERS 


Gun Plasma Density See Oeo mm 
Plasma Drift Velocity 7.0x10° m/s 
Delay time 1.5x10°° sg 
Time of flight 1.36x107° s 
Cathode radius 2.5x107° m 
Anode radius SAO 0s— am 
Gun-cathode distance 9.5x10°° m 
Switch length Al scOe an 


A. IONIZATION RATIO 
The switch parameters listed in Table 2 which will be 
varied during the evaluation are the density of the plasma 


guns and the delay time between the firing of the plasma guns 
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and the firing of the generator. In addition, the effect of 
the electron energy, monolayers desorbed and the density 
profile of the desorbed particles will be considered. 

I Density profile 

The density profile utilized is explained in 
Appendix A. Essentially, the expanding carbon cloud is 
divided into 5 concentric annular regions of equal spacing. 
Each region is assigned an equal number of neutral carbon 
particles. This method produces a slightly non-uniform 
density profile due to the cylindrical geometry. However, 
the profile can be altered by applying a weighting factor f. 
For f = 1, the number of particles in each region is equal. 
For 0.5 < £ < 1, the number of particles is more heavily 
weighted to the outer regions while 1 < f < 2, more particles 
are allocated within the inner regions. 

The effect of this particle weighting can be seen 
in Figures 2, 3, 4 and 5. For 50 eV electrons, the 
1onization of the carbon particles saturates within 40-45 
nanoseconds. Allocating more particles in the outer regions 
produces a less varied plasma density difference between 
regions than allocating them to the inner regions. By 
weighting the inner region, an order of magnitude plasma 
density difference can be achieved. 

Ze Desorbed layers 
As seen in Figure 6, variations in the number or 


fractional number of monolayers desorbed affects the 
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ionization ratio by the same factor. Increase the number of 
monolayers desorbed by a factor of 10, and a factor of 10 
increase in the ionization ratio is observed. This 
information provides a baseline to consider for the number of 
monolayers that must be removed in order to produce any 
Significant effect on the electric field or the sheath width 
as proposed by Schwirzke. Equations 1 and 2 scale with the 
square root of the plasma density. Figure 6 clearly 
demonstrates that in order to influence the electric field, 
at least 0.1 monolayers of carbon particles must be desorbed 
from the surface. In order to see the trends more clearly, 
one monolayer will be utilized for all future comparisons. 
oye Rlectron energy variations 

From Table 1, it is clear that electrons with 
energies below 100 eV have larger ionization cross sections 
than those of higher energies. Figure 7 demonstrates that 
the lower energy electrons ionize the carbon layer much more 
rapidly than the higher energy electrons. Based on the known 
time scales of the plasma opening switch, it can also be 
inferred that the voltage across the sheath must remain 
relatively low ( 20° S £ = 1008) durimeGe@eme initial conducere 
phase of the current. This allows for rapid plasma density 
increase near the cathode surface. However, the electric 
field in Equation 2 is dependent on the density and the 
applied voltage. If the anode-cathode potential difference 


remains low, then the density must increase significantly to 
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provide the electric fields necessary to sustain electron 
emission. 
4. Delay time 
By varying the delay time between the firing of the 
plasma guns and the generator, the size and density of the 
desorbed carbon cloud are modified. This results in an 
overall decrease in the ability of the electrons to ionize 
the carbon atoms as the delay time increases. As stated 
previously, the effect of the neutral carbon density on the 
mean free path is significant. Figure 8 demonstrates this as 
the density ratio D decreases as the delay time increases. 
oe Plasma gun density 
As the plasma gun density varies, the ionization 
ratio also varies. Increased plasma gun densities require 
that the ionization of the carbon atoms be much greater 
before the effect is noticeable. Whereas decreased plasma 
gun densities can produce a significant ionization ratio 


almost immediately. Figure 9 demonstrates this principle. 


Be PLASMA DENSITY GRADIENT EFFECT 

The plasma density gradient produces a pressure force in 
the vicinity of the cathode surface. The influence of the 
density gradient on the plasma particle motion is more 
difficult to assess with the given information. The data 
presented is static while the effect of the density gradient 


is dynamic in nature. 
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The carbon atoms forming ehem@meteral Cloudede nce ree 
the effect of the electric field. Therefore, the ionization 
of the neutral particles is not affected by the density 
gradient. The ionized particles, on the other hand, will be 
accelerated by the electric field toward the cathode or anode 
depending on the charge and direction of the field at that 
point in space. 

The electric field profile due to the density gradient 
will be of the form presented in Figure 10 where negative E 
represents the direction toward the cathode and positive E is 
toward the anode. The data for the electric field produced 
by the plasma density gradient is presented only between the 


cathode surface and the first annular region. 
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neeres Il i218, 14, and ls show that the electric 
field due to the density gradient rapidly saturates. The 
reason for this is that Equation 15 becomes constant when the 
ratio D becomes large for a fixed electron energy. Equation 
15 then becomes a function of the distance over which the 
gradient is calculated. This is essentially a function of 
the expansion time of the neutral carbon cloud as shown in 
Figure 14. In addition, the energy of the plasma also affect 


Biewelectric tileld of Equation 15 as shown in Figure 13. 
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Ve NCLUSIONS 


The purpose of this study was to conduct some 
preliminary testing of Schwirzke's Plasma Sheath Model in an 
effort to determine the effect of the neutral carbon 
ionization process on the electric field experienced at the 
cathode. Inherent in that attempt was the determination that 
ionization was indeed occurring. In addition, the effect of 
the plasma density gradient was considered. 

It was shown that even if modest levels of surface 
contaminants were desorbed, ionization was appreciable. Asa 
result of this study, a baseline level of desorbed particles 
was established. Although the results are approximations 
based on ideal conditions, this information may be helpful 
for future investigations. 

It was shown that variations in the delay time, plasma 
density and the desorption level provided the most dramatic 
effects on the density ratio D. This, in turn, affects the 
electric field and the electron field emission current drawn 
from the cathode surface. 

The most interesting effect of the neutral particle 
ionization, is the effect due to the plasma density gradient 
which was shown to produce significant electric fields near 
the cathode. This effect could prove to be significant in 


future studies for two reasons. First, the plasma density 
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fedeeeme Produces electric fields on either side of the 
1onizing region which are diametrically opposed. Second, the 
expansion of the ions is confined on one side by the cathode 
and essentially free on the anode side. 

The electric field produced by the density gradient near 
the cathode accelerates the ions toward the cathode while the 
electrons are essentially trapped within the opposing 
electric fields near the location of E = 0. As the ions are 
lost to the cathode surface, the electrons form a negative 
space charge layer which reduces the emission of electrons 
from the cathode. This could allow the current to be 
transfered to the load. At the same time, the anode side of 
the density gradient produces an electric field which opposes 
the electron motion toward the anode. Depending on the rate 
of expansion and ionization, this electric field could be 
responsible for the observed flow of ions near the anode. As 
the dense plasma cloud continues to expand, the electric 
field is reduced and the density ratio approaches unity 
again. 

The effects of the anode pressure gradient and rate of 
expansion were not included in this study. However, the 
dynamic effects of the density gradient and resultant 


electric fields could prove important. 
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APPENDIX A NBUTRAL DENSITY PROFILE 


The basic notion of the initial number and number 
density in each of the five concentric anular rings is 
presented in the body of Chapter II and Equations 5 and 6. 
The weighting factor is based on a straight line curve 
approximation where the number of neutral particles remains 
COnstant:, 

The inner most region is multiplied by the weighting 
factor £, which ranges as 0lo7s 22 serosa Al shewaee 
schematic of the annular sub-regions around the cathode. The 
equation of the line then depends on the slope and the y- 
intercept. The y-intercept is determined by the factor f 
multiplied by the number of particles in the innermost 
region. The slope is determined by forcing the middle region 
to always have the same number of particles as was originally 
allocated 

This method of weighting is somewhat restrictive 
althougth for the purposes of this study, the factor f is 
easily identified. For example, for f = 1, each region has 
the same number of particles while f = 2, immediately alerts 
the observer that the inner region is heavily weighted with 


Pareue les 
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APPENDIX B PROGRAM LISTING 


program sheath (input, output); 


{The purpose of this program is to calculate the ionization 
of neutral particles which are desorbed from a cathode 
surface in a plasma opening switch configuration. It also 
calculates the electric field produced by the density 
gradient. Calculations used are based on fundamental plasma 
physics equations where no attempt has been made to apply 
kinetic theory. } 


uses SANE; 


COnseE 
charge = 1.681, 
emass = 9.1E-31; 
imass = 2.004E-26; 
Permi %="35 coea te, 
cathradius = 2.5E-2; 


anoderadius = 5E-2; 
Grifievel = 7Ee, 
Evel 250 


Swwidth = 4E-2; 
swlength = 9.5E-2; 
delay = 1.5E-6; 
londensity = 8E19; 


var 
voltage, ivelocity, evelocity,current,selvolt,neutdensity, 
lamda, Ssigma,ncur,endcur, catharea, anodearea,eondensity, 
11, Je, xpcatlo,@eutradius, 
stepcur,irad,n, length,eflux,neutpart, 
partsinregion,f,left,slope,parts,lam, ionized:real; 
idens,deltat,dist,ri,ro,pos:extended; 
tyi,steps, kay mst ee eger; 
filel:text; 
Chie clidaar 
nparticles,iparticles,ndensity, 
idensity,ep:array [12.6 |Metsreau, 
match:boolean; 


procedure setcon; {set the initial constants for the given 
parameters } 
begin 

ivelocity:=sqrt (2*charge*voltage/imass) ; 
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evelocity:=sqrt (2*charge*voltage/emass) ; 
neutradius:=cathradius+ (delay-swlength/driftvel) *tvel; 
X:=neutradius-cathradius; 
@ast :=x/5; 
Geemarea:—Z2Z*p1i~cCacthnradius*swwidth; 
anodearea:=2*pi*anoderadius*swwidth; 
Houle paGhs—-caenarean*].Z217E19; 

end; 


procedure cathodedensity; {determine the current density and 
electron flux} 


begin 
jJi:=londensity*charge*driftvel; 
je:=(current- (ji*catharea) ) / (anodearea) ; 


eflux:=je/charge*1lE-9;if eflux<=0 then eflux:=1; 
{part/square meter nanosecond} 
end; 


procedure getinfo; {assume a straight line current profile} 
begin 
write('Enter the anode-cathode voltage '); 
readin (voltage) ; 
write('Enter the number of monolayers to remove '); 
mean (nh) ; 
Wace ('’EmtLer the starting current in kiloamps "); 
mean (mCcur) ; 
Gi@ment > =—neur*i1b3; 
write('Enter the ending current in kiloamps '); 
readin (endcur) ; 
endcur:=endcur*1E3; 
write ('Enter the number of steps to take '); 
readin (steps); 
stepcur:=(endcur-current) / (float (steps) ); 
end; 


procedure setsigma; {assign the proper cross section} 
var 
setvolt:integer; 
begin 
setvolt:=Round (selvolt); 
mecch;=DRUE: 
case setvolt of 
Z0-wSigma:—1b-20; 
SOs siama:—1.8E—-20; 
a0) sicma:—=2.1E-20; 
a0) Sigma :—=Z35520 > 
60: sigma:=2.4E-20; 
moe sroame:—2.3E-20; 
Stee stagma:=Z.2E-20; 
See stoma:=2Z.15E-20; 
HOG Saamas—=Z.1E-20; 
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200: Siomagalsh—20-; 
300: sigmas—1. 2n—20- 
400: sigma:=9E-21; 
500: sigma:=8E-21; 
600: sigma:=7E-21; 
700: sSigma:=6-5E7Z1i- 
800: sigma:=6E-21; 
S10 5 Sse ys) 5 Sie AIS 
1000: Sigma:=Se-21; 

otherwise 
match:=FALSE; 

end; 

end; 


procedure profile; {write the equation of the line} 
begin 
if f<>l then begin 
left:=f*partsinregion; 
slope :=2* (1=£) *partsinregiony =; 
end; 
if f=1 then begin 
left:=partsinregion; 
slope:=0; 
end; 
end; 


procedure openfiles; 
begin 

rewrite (filel, 'data set 1'); 
end; 


begin 

openfiles; 

match:=FALSE; 

while not match do begin 
getinfo; 
selvolt:=voltage; 
setsigma; 
end; 


setcon;length:=dist; 

partsinregion:=neutpart/10; 

writeln('Enter the factor ');readln(£)7i1£ f-0¢emen halt; 
DrOfELIe) pos sa Gais ey 


write (filel, 'voltage = ',voltage:11); 
write (ftilel,’dlengeh =) 'Fx- ll) 
write(filel, 'monolayer = ',n:11); 
Wrelteelnd(tilel Lmeactor =m eee 
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for k:=1l to 5 do begin {assign the particles to the regions} 
nparticles[k] :=slope* (pos-dist/2)+left; 
iparticles[k] :=0; 
pos:=post+dist; 
Bice! 


tetas —1) to steps do begin 
Cacthnodedensity; rl:=cathradius; ro:=cathradius+dist; 


for j:=1 to 5 do begin 

neéenstty |) |) :—mparticles ([j]// (p1*swwidth* (sqr (ro) - 
Se (tei) ) ); 

lamda:=1/ (ndensity[ 4] *sigma) ; 

ionized:=(current/charge-ji) *1E-9* (1-exp ( (- 
length) /lamda) ) ; 

Llparticles[j] :=iparticles[ j]+ionized; 

tego levid| :—-lLparticles([j)/ (pi*swwidth* (sqr(ro)- 
Sciae!)°) )\ 7 

nparticles[j]:=nparticles[jJ-iparticles[j]; 

if nparticles[{j]<=0 then nparticles[j] :=1; 

Pes —Gere eo: —=rotdist; 


end; 
if 1 mod 5 = 0 then begin {print out every 5 nanoseconds} 
write (filel,i:5);write(filel,' '); 


write (filel,voltage:11); 
Wiese (fi lel, curréent:11) > 
Poel —ieco 5 do begin 
write (filel, (idensity[m]+iondensity) /iondensity:11); 
end; 
Vieni ohr (9) ) ? 
ep[1]:=abs((2*selvolt) /(2*iondensity+idensity([1])* 
Gidenouse,| 1) /dist)> 
write (filel,ep[1]:11); 
home Mo=2et0 5 do begin 
ep[m] :=abs ((2*selvolt) / (idensity[ml]+idensity{[m]+ 
2 Pemeemorey)) ~(1density({m)]—-idensity[m-1]) /dist); 
menee(fehe ls, ep[(m] 2:11) > 
end; 


ep[6] :=abs ((2*selvolt) / (idensity[5]+2*iondensity) *idensity[5] 
Foams ):; 
Wepeclm(filel,ep[6] 2:11) ; 
end; 


GUrrent:—currentt+stepcur; 
end; 

ebose(tilel) s 

readin; 

end. 
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